RAISHIN: A High-Resolution Three-Dimensional General Relativistic Magnetohydrodynamics Code

We have developed a new three-dimensional general relativistic magnetohydrodynamic (GRMHD) code, RAISHIN, using a conservative, high resolution shock-capturing scheme. The numerical fluxes are calculated using the Harten, Lax, & van Leer (HLL) approximate Riemann solver scheme. The fluxinterpolated, constrained transport scheme is used to maintain a divergence-free magnetic field. In order to examine the numerical accuracy and the numerical efficiency, the code uses four different reconstruction methods: piecewise linear methods with Minmod and MC slope-limiter function, convex essentially non-oscillatory (CENO) method, and piecewise parabolic method (PPM) using multistep TVD Runge-Kutta time advance methods with second and third-order time accuracy. We describe code performance on an extensive set of test problems in both special and general relativity. Our new GRMHD code has proven to be accurate in second order and has successfully passed with all tests performed, including highly relativistic and magnetized cases in both special and general relativity. Subject headings: accretion, accretion disks black hole physics magnetohydrodynamics: (MHD) method: numerical -relativity National Space Science and Technology Center, 320 Sparkman Drive, VP 62, Huntsville, AL 35805, USA; [email protected] Center for Space Plasma and Aeronomic Research, University of Alabama in Huntsville Department of Physics, Kumamoto University, Kurokami, Kumamoto, 860-8555, Japan Department of Physics and Astronomy, The University of Alabama, Tuscaloosa, AL 35487, USA NASA-Marshall Space Flight Center, National Space Science and Technology Center, 320 Sparkman Drive, VP 62,Huntsville, AL 35805, USA NASA Postdoctoral Program Fellow/ NASA Marshall Space Flight Center